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Abstract We review the various techniques through which wind properties of mas- 
sive stars - O stars, AB supergiants, Luminous Blue Variables (LBVs), 
Wolf-Rayet (WR) stars and cool supergiants - are derived. The wind 
momentum-luminosity relation (e.g. Kudritzki et al. 1999) provides a 
method of predicting mass-loss rates of O stars and blue supergiants 
which is superior to previous parameterizations. Assuming the theo- 
retical Z 0,5 metallicity dependence, Magellanic Cloud O star mass-loss 
rates are typically matched to within a factor of two for various cali- 
brations. Stellar winds from LBVs are typically denser and slower than 
equivalent B supergiants, with exceptional mass-loss rates during giant 
eruptions M = 1(T 3 . . . W' 1 M Q yr' 1 (Drissen et al. 2001). Recent 
mass-loss rates for Galactic WR stars indicate a downward revision of 
2-4 relative to previous calibrations due to clumping (e.g. Schmutz 
1997), although evidence for a metallicity dependence remains incon- 
clusive (Crowther 2000). Mass-loss properties of luminous (> 10 5 Lq) 
yellow and red supergiants from alternative techniques remain highly 
contradictory Recent Galactic and LMC results for RSG reveal a large 
scatter such that typical mass-loss rates lie in the range 10 -6 . . . 10~ 4 M e 
yr _1 , with a few cases exhibiting ~ 10 _3 Mq yr _1 . 

1. INTRODUCTION 

Winds are ubiquitous in massive stars, although the physical pro- 
cesses involved depend upon the location of the star within the H-R 
diagram. Mass-loss crucially affects the evolution and fate of a massive 
star (Meynet et al. 1994), while the momentum and energy expelled 
contribute to the dynamics and energetics of the ISM. The interested 
reader is referred to Lamers & Cassinelli (1999) for the topic of stellar 
winds in general, or Kudritzki & Puis (2000) for a detailed discussion 
of mass-loss from OB stars. This review will cover the broad topics of 
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mass-loss from early-type (OBA, Wolf-Rayet, Luminous Blue Variable) 
and late- type (yellow/red supergiants) stars. Theoretical predictions for 
hot star winds follow line-driven radiation pressure (Castor et al. 1975; 
Pauldrach et al. 1986). Pulsations are thought to initiate the winds of 
cool supergiants, together with other mechanisms (e.g. Alfven waves). 
Radiation pressure on dust grains helps to maintain outflows in their 
cool, outer envelopes (Wickramasinghe et al. 1966). 

2. HOT STAR WIND DIAGNOSTICS 

Global wind properties can be characterized by mass-loss and wind 
velocity. In hot stars, the former depends on application of varying 
complexity of theoretical interpretation, while the latter can fortunately 
be directly measured with minimal interpretation. 
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Figure 1 Comparison between HST UV and FUSE far-UV spectra of Sk80 (07Iaf) 
with synthetic spectra (including shocks) at T e g= 37,500K and 32,000K, illustrating 
the improved agreement for the latter value (Fullerton et al. 2000) . 
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2.1. UV AND FAR-UV SPECTROSCOPY 

Ultraviolet P Cygni profiles, ubiquitous in O-type stars (Walborn et 
al. 1985) provide a direct indication of stellar winds. Such metal reso- 
nance lines are generally analysed via the Sobolev with exact integration 
(SEI) method (Lamers et al. 1987; Haser et al. 1998). Mass-loss de- 
terminations require knowledge of elemental abundances, the degree of 
ionization of the ion producing the line, and an assumed form of the 
velocity law that is predicted by radiation driven winds. In practice, 
accurate wind velocities can be readily obtained from saturated UV P 
Cygni profiles (Groenwegen et al. 1989; Prinja et al. 1990), while unsat- 
urated P Cygni profiles provide a better means of determining mass-loss 
rates, albeit still subject to ionization results from non-LTE atmospheric 
models. Such models are subject to uncertain contributions of ionizing 
radiation from shocks in the stellar outflow. In view of these difficul- 
ties, Lamers et al. (1999) have combined column densities from unsat- 
urated UV P Cygni profiles with independently determined mass-loss 
rates (from radio or Ha, see below) to derive empirical ionization frac- 
tions. 

The Far-Ultraviolet Spectroscopic Explorer (FUSE) offers the possi- 
bility of determining mass-loss rates solely from UV data alone, since 
the far-UV significantly increases the number of stellar wind diagnos- 
tics, including Cm, Niii-iv, Ovi, Siv-vi and Pv. The additional 
ionization information provided by FUSE permits stellar temperature 
determinations. Fullerton et al. derived T e g=32-33kK for two Magel- 
lanic Cloud 07 supergiants from UV wind metal diagnostics, in contrast 
to T e g=38kK from the usual optical helium plane-parallel hydrostatic 
methods (see Fig. 1). Although previously derived mass- loss rates were 
confirmed, lower bolometric corrections implied reduced stellar luminosi- 
ties, with consequences for current mass-loss calibrations. 

2.2. OPTICAL AND NEAR-IR 
SPECTROSCOPY 

Although Ha has long been recognized as the prime source of mass- 
loss information in early-type stars, accurate determinations rely upon 
a complex treatment of the lines and continua, i.e. non-LTE, spherically 
extended models which treat the sub- and supersonic atmospheric struc- 
ture. Use of such codes also provide surface temperatures and gravity in 
OB stars via optical photospheric lines (e.g. Hei-n in O stars). Prob- 
lems with using fits to Ha to derive mass-loss rates include blending 
with Hell A6560 in O stars, together with nebular contamination from 
Hn regions (e.g. NGC346; Walborn et al. 1995) and broadening of the 
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central emission component by stellar rotation. In A supergiants, Ha 
behaves as a scattering line, with a characteristic P Cygni profile, so 
they have the advantage over other hot stars that wind velocities and 
mass-loss rates may be simultaneously determined in external galaxies 
without the need for UV spectroscopy 

Although moderate to high spectral resolution observations of OB 
stars are rare at near-IR wavelengths (e.g. Fullerton &; Najarro 1998), 
analogous lines to optical wind (and photospheric) diagnostics are avail- 
able (Bohannan & Crowther 1999). This provides the prospect of the 
determination of mass-loss properties for hot stars obscured at optical 
wavelengths via IR diagnostics (Najarro et al. 1994). 

2.3. IR-RADIO CONTINUA 

Winds in hot stars can be readily observed at IR-mm-radio wave- 
lengths via the free-free (Bremsstrahlung) 'thermal' excess caused by 
the stellar wind, under the assumption of homogeneity and spherical 
symmetry. Mass-loss rates can be determined via application of rela- 
tively simple analytical relations (Wright & Barlow 1975; Panagia & 
Felli 1975). Barlow & Cohen (1977) have used IR excess fluxes to de- 
termine mass-loss rates for a large sample of Galactic OBA supergiants, 
which generally compare favourably with recent Ha results for A and 
B supergiants (Kudritzki et al. 1999). Unfortunately, OB stars with 
relatively weak winds do not show a strong IR excess or radio flux, so 
mass-loss results have solely been for nearby hot stars with dense winds 
(Bieging et al. 1989; Drake Sz Linsky 1989). Observations at multiple 
frequencies are necessary to ensure against non-thermal (synchrotron) 
radio emission from colliding winds (e.g. Chapman et al. 1999). Ha 
and radio continuum mass-loss rates agree reasonably well (Lamers Sz 
Leitherer 1993) despite sampling quite different parts of the stellar wind 
- ~1.5 stellar radii for Ha and ^1000 stellar radii for the radio photo- 
sphere - placing constraints on relative clumping factors in these regions. 

3. RESULTS FOR OBA STARS - ROLE OF 
METALLICITY? 

Although the techniques used to study hot, massive stars are similar, 
we first discuss 'normal' early-type stars, and leave 'exotica' until subse- 
quent sections. The mass- loss rate of a star is a fundamental ingredient 
in evolutionary models. Unfortunately, radiatively driven wind theory 
in its current form (Puis et al. 1996; Vink et al. 2000) is unable to 
reproduce the observed mass-loss properties for hot luminous stars suffi- 
ciently accurately (see below). Nevertheless, since their winds are driven 
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by photon momentum transfer through metal line absorption, the stel- 
lar momenta and wind velocities are expected to be functions of stellar 
metallicity, Z, such that mass-loss rates scale with Z 0,5 for O stars and 
early B stars, Z ' 8 for mid-B supergiants, and Z 1 ' 7 for A supergiants 
(Kudritzki k Puis 2000). 

3.1. TERMINAL WIND VELOCITIES 

Large compilations of wind velocities, Voo, of O-stars and AB super- 
giants have been made by Prinja et al. (1990) and Lamers et al. (1995), 
revealing a spectral type and luminosity class dependence. Lamers 
et al. highlighted the so-called 'bi-stability' jump in wind properties 
around B1.5 (~21,000), above which ~ 2.65-u e sc (escape velocity), 
and -Woo ~ 1.4i; csc below, resulting from the change in ionization of the 
elements contributing to the line force (Vink et al. 1999). Wind veloci- 
ties of LMC O stars differ little from Galactic counterparts (Garmany & 
Conti 1985), while a more prominent effect is observed in the SMC, par- 
ticularly amongst early O stars (Walborn et al. 1995; Prinja k, Crowther 
1998). 

Table 1 Coefficients for wind momentum iuminosity reiationship of Gaiactic O-stars 
and AB supergiants (Kudritzki & Puis 2000) 



Sp Type 



log D 



x 



O dwarf /giant 
O super giant 
B0-1 supergiant 
B 1.5-3 supergiant 
AO-3 supergiant 



19.87±1.21 
20.69±1.04 
21.24±1.38 
17.07±1.05 
14.22±2.41 



1.57±0.21 
1.51±0.18 
1.34±0.25 
1.95±0.20 
2.64±0.47 



3.2. MASS-LOSS RATES 

Observationally, estimates of mass- loss rates from Ha, UV wind lines 
or radio fluxes have been available for the past 20 years, although is 
only recently that a reasonable degree of consistency has been reached 
between these different approaches via theoretical improvements. Puis 
et al. (1996) and Kudritzki et al. (1999) identified that O-stars and AB 
supergiants obey a wind momentum-luminosity relationship, as follows 

logMv^R/Ref 5 = logD + x\og(L/L Q ) 
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Stars of different spectral types have differing Dq and x values, listed 
in Table 1, because of the change of lines driving the stellar wind with 
T e s- 

Hot stars in the Magellanic Clouds provide us with the means to 
compare emirical mass-loss properties with predictions at lower Z (e.g. 
Garmany & Conti 1985). This has proved to be difficult given (i) the 
inability to measure radio fluxes for extra-galactic hot stars; (ii) the 
need to combine high quality UV and optical spectroscopy with sophis- 
ticated model techniques; (iii) the limited metallicity range spanned, 
0.15-1^0. Although the number of stars analysed in detail for the Mag- 
ellanic Clouds remains small (Puis et al. 1996), the mass- loss rates of 
LMC and SMC O stars tend to be lower than Galactic counterparts. As 
an example, let us consider Sk80 (07Iaf+) in the SMC giant Hn region 
NGC 346 which has a metallicity of 0.2Z Q (Haser et al. 1998; Peim- 
bert et al. 2000). From the wind momentum-luminosity relationship we 
would expect an equivalent Galactic supergiant to have a mass loss rate 
four times higher than its observed value (Puis et al. 1996), suggesting 
a higher exponent of Z°' S in this case. 

Table 2 Comparison (M(T h, L)/M b s ; cr in parenthesis) between predicted OBA 
mass-loss rates from various parameterizations with empirical data from Puis et al. 
(1996) and Kudritzki et al. (1999). Parameterizations are de Jager et al. 1988; 
Lamers & Cassinelli 1996; Vanbeveren et al. 1998; Kudritzki & Puis (2000); Vink 
et al. (2000) or Achmad et al. (1997) for A supergiants. Except for Lamers & 
Cassinelli, no metallicity dependence is accounted for, so we here assume a metallicity 
dependence of Z - 5 for O stars, with 0AZ e for the LMC, 0.15Z Q for the SMC. 



Sped. 


Galaxy de Jager Lamers & Vanbeveren Kudritzki 


Vink 




type (#) 




et al. 


Cassinelli 


et al. 


k Puis 


et al. 




O (23) 


Gal 


0.9 (0.7) 


1.1 (0.6) 


0.9 (0.7) 


1.3 (0.8) 


1.8 


(1.3) 


O (6) 


LMC 


0.5 (0.5) 


1.3 (0.7) 


0.5 (0.3) 


0.7 (0.3) 


1.3 


(0.4) 


O (5) 


SMC 


1.3 (0.9) 


2.5 (1.5) 


1.3 (0.9) 


2.3 (1.3) 


6.5 


(7.9) 


B(14) 


Gal 


6.8 (4.6) 


11.3 (7.7) 




1.0 (0.3) 


38.8 (29.4) 


A (4) 


Gal 


5.0 (1.6) 


3.7 (1.9) 




1.0 (0.3) 


2.0 


(1.0) 



Evolutionary calculations require M(T e g,L) parameterizations based 
on empirical data. Up until recently, the only compilation covering hot, 
luminous stars was that of de Jager et al. (1988). In an attempt to up- 
date the mass-loss calibration for evolutionary calculations, Vanbeveren 
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et al. (1998) provided an updated relationship for Galactic OB stars, 
while Lamers & Cassinelli (1996) additionally considered LMC/SMC O 
stars. Most recently, Achmad et al. (1997) and Vink et al. (2000) have 
provided a set of predictions from current radiative driven winds for A, 
F and G supergiants and OB stars, respectively. How do these different 
parameterizations compare with observed mass-loss rates? 

We have taken empirical data, M Q b s , from Puis et al. (1996) for O 
stars and Kudritzki et al. (1999) for AB supergiants. In Table 2 we 
compare the mean ratio M(T c r , L)/M G b s for each calibration. Perhaps 
surprisingly, de Jager et al. (1988) fares as well as the recent Lamers 
& Cassinelli or Vanbeveren et al. (1998) parameterizations for O stars. 
In contrast with previous theoretical results, Vink et al. (2000) predicts 
sufficiently strong winds, albeit a factor of 2 too strong for Galactic 
O stars. Overall, the massdoss rates of AB supergiants are typically 
overestimated (by 5-10) with the exception of Kudritzki et al. (1999) 
and Achmad et al. (1997). Predictions by Vink et al. (2000) show the 
poorest agreement for mid-B supergiants. 

Overall, solely the wind momentumduminosity prescription provides 
excellent agreement with all OB A observations. One note of caution is 
warranted, since the T c g-calibration for O stars, and hence the luminos- 
ity scale, relies almost universally on results from plane-parallel studies 
which neglect stellar winds (recall Sect 2.1). 

3.3. ROTATION 

O star rotational velocities (i>sini) are typically 100 kms -1 (Howarth 
et al. 1997), with important consequences for their evolution (Meynet 
& Maeder 2000). Friend & Abbott (1986) first considered the effect 
of centrifugal acceleration on hot star winds, extended by Bjorkman 
& Cassinelli (1993) to allow for the azimuthal dependence, who found 
that M increases (and Voo decreases) towards the equator, leading to an 
oblate 'wind-compressed disk'. 

Subsequently, Owocki et al. (1996) established that quite different 
predictions may result, by additionally allowing for non-radial compo- 
nents of the line force, and the increased polar radiation flux via gravity 
darkening (von Zeipel 1924). Due to a reduced escape velocity, the equa- 
torial velocity law is slower than at the pole, inhibiting the equatorial 
disk and the original scaling of mass-loss can be reversed, such that a 
prolate geometry results. As an example, a mid-B dwarf rotating at 
85% of its critical velocity was found by Petrenz & Puis (2000) to have 
a strongly prolate structure (p po le/Peq <15) with a polar /equatorial ter- 
minal velocity of 1030/730 kms -1 . The global mass-loss rate was not 
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found to deviate from its ID value by greater than 10-20%, except for 
supergiants close to the Eddington limit. 

3.4. STRUCTURE 

Observationally, evidence for structure in O stars winds has taken 
various forms, including absorption variability observed from optical 
(Fullerton et al. 1996) and UV (Prinja 1992) spectroscopic intensive 
monitoring, and the presence of soft X-ray emission (Chlebowski et al. 
1989), attributed to shocks in their winds (Lucy & White 1980). This 
view has recently been questioned by Chandra observations of £ Ori 
(Waldron Sz Cassinelli 2000) which reveal that the X-ray emitting gas is 
exceptionally dense, with no evidence for expansion. 

Theoretical studies of hot star winds reveal that line driven accelera- 
tion is subject to a strong instability, causing a variable and structured 
winds (Owocki et al. 1988; Feldmeier 1995). Since the region where 
Ha forms in OB stars (< 1.5i?*) is not anticipated to be heavily struc- 
tured, mass-loss determinations may not be greatly affected by clump- 
ing/structure. Recent extensive monitoring of HD 152408 (WN9ha or 
08Iafpe) by Prinja et al. (2000) indicate variations of ±10% in global 
mass-loss rate, although with a time averaged Ha profile that is remark- 
ably stable over several years. 

4. LUMINOUS BLUE VARIABLES 

The stellar wind properties of LBVs are poorly constrained, although 
it is apparent that their winds are generally slower and denser than OBA 
supergiants (Crowther 1997), with little variation in mass- loss during 
their 'normal' irregular excursions across the H-R diagram (Leitherer 
1997). To illustrate this, let us consider P Cygni, whose parameters 
have been determined by Najarro et al. (1997). Compared to normal 
B1.5 supergiants (780 kms -1 ; Lamers et al. 1995), P Cyg has a very slow 
wind (185 kms -1 ) and a mass-loss rate which exceeds the Galactic wind 
momentum relation (Sect. 2) by a factor of 15. Recall that from Petrenz 
& Puis (2000), non-rotating models of supergiants close to the Eddington 
limit may overestimate global mass- loss rates by a factor of ~2, so that 
differences may be somewhat less severe. There are numerous examples 
of extreme supergiants which are not LBVs yet show remarkably high 
mass- loss rates and low wind velocities (e.g. Hillier et al. 1998). 

In the case of giant LBV eruptions, information is very scarce and 
up until recently solely reliant on the 1840's rj Car eruption causing 
the formation of the Homunculus, in which ~3 M was ejected over 
~30 years (Davidson & Humphreys 1997), corresponding to 0.1 M Q 
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yr . Within the past few years, two extra-galactic LBVs have had 
sudden giant eruptions, HD 5980 in the SMC (Barba et al. 1995) and 
VI in NGC2363 (Drissen et al. 1997). These eruptions, in stars at low 
metallicity environments, appear to be quite different from that of r] Car 
during its eruption, since peak mass-loss rates of each are ~ lO -3 M yr _1 
(Drissen et al. 2001), equivalent to the current mass-loss for i] Car 
(Hillier et al. 2000). 

5. WOLF-RAYET STARS 

The diagnostics of mass-loss for Wolf-Rayet stars mimic O stars ex- 
cept that their stronger winds mean that UV/optical wind diagnostics 
are not restricted to resonance lines or Ha. Wind velocities can be accu- 
rately determined either from the usual UV P Cygni profiles (Prinja et 
al. 1990), or from optical/IR Hei P Cygni profiles (Eenens & Williams 
1994). Radio and mm fluxes have been used to determine mass-loss 
rates of Galactic WR stars (Leitherer et al. 1997), while recent theoreti- 
cal progress (Schmutz 1997; Hillier & Miller 1998) now permits detailed 
non-LTE analysis of UV and optical observations. WR winds are typi- 
cally spherical (Harries et al. 1998), but clumped (Lepine et al. 2000). 
Spectroscopic tools make use of electron scattering wings in WR winds 
to estimate clumping factors since their strength depends linearly on 
density, versus the square of density for the underlying recombination 
line (Hillier 1991). Typical volume filling factors are 5-10%, such that 
mass-loss rates in a clumpy medium are 3-4 times lower than those from 
homogeneous models (e.g. Hillier & Miller 1999). 

Although great progress has been made in recent years, detailed analy- 
sis requires labour intensive use of complex models, so the sample of stars 
studied in detail remains small. Recently, Nugis & Lamers (2000) com- 
bined empirical IR-radio observations with analytical models accounting 
for clumping to estimate mass-loss rates as a function of luminosity, He 
(Y) and metal (Z) mass fractions, i.e. 

log(MwN) = -13.6 + 1.63 log L/L Q + 2.22 log Y 

log(M W c) = -8.3 + 0.84 log L/L Q + 2.04 logY + 1.04 log Z 

Recent evolutionary calculations for hydrogen-poor WR stars follow 
mass- loss calibrations of Langer (1989) or Vanbeveren et al. (1998). 
Table 3 compares predictions from these calibrations with results from 
detailed non-LTE line blanketed spectroscopic analyses, assuming the 
Schaerer & Maeder (1992) mass-luminosity relationship. The formula- 
tion of Langer (1989) overestimates mass- loss rates by factors of 2-4, 
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while Vanbeveren et al. (1998) and especially Nugis & Lamers (2000) 
show much improved agreement with spectroscopic results. 

Table 3 Comparison between WR mass-loss rates, in 10 _5 MQyr _1 , predicted by 
Langer (1989, L89), Vanbeveren et al. (1998, VDV98), semi-empirical calibration of 
Nugis & Lamers (2000, NL2000) with detailed non-LTE spectroscopic results (nLTE) 
of hydrogen poor Galactic WN and WC stars. 



Star Sp Type L89 VDV98 NL2000 nLTE Ref 



WR6 WN4 


12.5 


5.5 


5.7 


3.2 


a 


WR147 WN8(h) 


9.5 


4.5 


3.2 


2.5 


b 


WR111 WC5 


2.8 


2.0 


1.3 


1.5 


c 


WR90 WC7 


4.8 


3.2 


2.6 


2.5 


d 


WR135 WC8 


2.0 


1.6 


1.6 


1.2 


d 



(a) Schmutz (1997); (ft) Morris et al. (2000); (c) Hillier & Miller (1999); (d) Dessart et al. 
(2000) 

The question of a Z-dependence of mass-loss in WR stars remains 
open. There is a clear distinction amongst spectral subtypes in galaxies 
of differing metal content, which has been qualitatively explained from 
evolutionary and spectroscopic models (Smith & Maeder 1991; Crowther 
2000). LMC WN stars (Crowther k Smith 1997; Hamann & Koesterke 
2000) show negligible spectroscopic difference from Galactic counter- 
parts. The situation is less clear for the SMC, since most WR stars are 
complicated by binarity. Sk41 is a possible exception, whose properties 
compare closely with Galactic or LMC counterparts (Crowther 2000). 
Further progress requires detailed study of apparently single WR stars 
beyond the Magellanic Clouds (e.g. Smartt et al. 2000), especially in 
nearby metal poor galaxies such as IC10 (Massey k Armandroff 1995). 

6. MASS-LOSS DIAGNOSTICS IN YELLOW 
AND RED SUPERGIANTS 

Various techniques have been applied to quantify mass-loss in yellow 
(FG) and red (KM) supergiants, whose wind driving mechanism remains 
uncertain. UV/optical observations of binaries with an early- type dwarf 
companion permit analysis of the red-supergiant chromosphere and wind 
in absorption (Bennett et al. 1996). Radio continua measurements have 
provided upper limits to mass-loss rates of F-supergiants (Drake k Lin- 
sky 1986) and M-supergiants (e.g. Harper et al. 2001). Alternatively, 
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sub-mm CO emission or mid-IR dust emission permit mass-loss deter- 
minations in cool supergiants. CO emission is a reliable probe for AGB 
stars, although the UV radiation field from supergiant chromospheres 
could lead to a different CO emission per unit mass loss from star to 
star. Radial pulsations are generally assumed to provide the mechanism 
of initiating mass- loss in M supergiants, as evidenced from circumstel- 
lar dust shells (Bowers et al. 1983), with some notable exceptions - 
VY CMa (M5Ia) and p Cas (F8Ia+). Outflow velocities of RSG are 
typically 10-40 kms -1 (Jura & Kleinmann 1990). 
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Figure 2 Mass-loss rates derived for 
Galactic (open circles) RSG by Jura & 
Kleinmann (1999) and LMC (filled cir- 
cles) RSG by Reid et al. (1990) us- 
ing the Jura (1986) formula, including 
predictions from de Jager et al. (1988, 
solid) and Reimers (1975, dot-dash) for 
Tgs =3,500K. Fits to Reid et al. data are 
by Vanbeveren et al. (1998, dashed) and 
Salashnich et al. (1999, dotted). 



— <^ 


1 ' 








• 


-3 


IRAS04530 
• 


~ 7104 IRAS04553-6825 / 


-4 


" 8 


- — / ° 


-5 


• 


• • 




/ • 

so 




-6 


a Ori 


Gal (Sylvester) o - 
LMC (van Loon) • 




i i 






5.0 5.2 5.4 5.6 5.8 
log L/L 

Figure 3 Mass-loss rates derived from 
modelling mid-IR spectroscopy of Galac- 
tic (open circles) and LMC (filled circles) 
RSG by Sylvester et al. (1998) and van 
Loon et al. (1999). Distances to Galac- 
tic RSG follow Sylvester et al. except for 
W Per and p Cep, which have been ad- 
justed to distances derived by Garmany 
& Stencel (1992). Lines are as described 
in the adjacent caption. 



Since we are concerned with massive stars in this review, we shall re- 
strict our discussion to the most luminous 'hypergiants' (de Jager 1998), 
with log (L/Lq) >5. Studies have largely focussed on M supergiants, 
with only highly unusual FG supergiants studied in detail (e.g. p Cas, 
HR8752, IRC+10420). The latter is widely interpreted as a RSG rapidly 
moving blueward (Jones et al. 1993), which has recently experienced 
huge variability in mass-loss rate, 10 -2 to 10~ 4 M Q yr _1 , via analysis of 
its circumstellar dust shell (Blocker et al. 1999). 
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Jura & Kleinmann (1990) applied the 'approximate' Jura (1986) for- 
mula for intermediate mass stars to IRAS 60/mi photometry of nearby 
yellow and red supergiants to estimate mass- loss rates. Reid et al. (1990) 
followed the same approach for a sample of LMC RSG, which have been 
used by Vanbeveren et al. (1998) and Salasnich et al. (1999) to pro- 
vide 'empirical' calibrations for evolutionary calculations, as illustrated 
in Fig. 2. Jura & Kleinmann (1990) data for Galactic RSG reveal sys- 
tematically lower rates, except for VY CMa and NML Cyg. Predictions 
from the empirical calibration of de Jager et al. (1988, solid) and the for- 
mulation of Reimers (1975, dot-dash) are also illustrated, for an assumed 
T e ff=3,500K (~M3I). Although the scatter is large, most Galactic data 
lie between the Reimers (1975) and de Jager et al. (1988) calibrations. 

Recently, van Loon et al. (1999) have modelled ISO spectroscopy and 
photometry of a similar LMC sample to Reid et al., and found rates 
which are 5-10 times lower. Meanwhile, Sylvester et al. (1998) have 
used the 10/xm silicate emission feature for a sample of Galactic M su- 
pergiants, following the Skinner & Whitmore (1988) relation, which was 
calibrated against sub-mm CO determinations. Their determinations 
revealed higher mass- loss rates than Jura & Kleinmann (1990). From 
this combined dataset, presented in Fig. 3, there is little evidence of 
a luminosity dependence of mass- loss. The Vanbeveren et al. (1998) 
and Salasnich et al. (1999) calibrations match Galactic data reasonably 
well, except for a Ori (M2.2Iab), which probably has the most secure 
determination. In contrast with Fig. 2, the majority of the LMC data in- 
dicate lower mass- loss rates than Galactic RSG, hinting at a metallicity 
dependence. However, van Loon et al. identified two LMC RSG with ex- 
ceptionally high mass-loss rates of order lO _3 M0yr _1 (IRAS 04553-6825 
and IRAS 04530-7104) which are unmatched in our Galaxy. 

Overall, the level of consistency between different IR/sub-mm tech- 
niques remains embarrassingly poor. Sylvester et al. (1998), using the 
power in W/im silicate emission (calibrated against earlier CO line data), 
derived a mass-loss rate for /i Cep (M2 la) which was 25 times higher 
than Jura & Kleinmann (1990) obtained from 60/im dust emission, or 
250 times higher than Josselin et al. (2000) based on an alternative CO 
line calibration of Olofsson et al. (1993). Indeed, VY CMa has values 
ranging over a factor of twenty - 3.5 x lO _4 M yr _1 (Stanek et al. 1995) 
to 1.6xlO" 5 M yr" 1 (Josselin et al. 2000). 

In summary, the mass- loss rates for RSG are by far the most contro- 
versial when recommending empirical calibrations for evolutionary cal- 
culations. The majority of mass-loss estimates lie between the Reimers 
(1975) and Vanbeveren et al. (1998) calibrations. This is illustrated 
in Table 4, where we compare rates assumed in the recent evolutionary 
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model of Meynet & Maeder (2000) for an initially rapidly rotating 40M Q 
star with 'empirical' calibrations of Kudritzki Sz Puis (2000) for OBA 
stars, Vanbeveren et al. (1998) for RSG and Nugis Sz Lamers (2000) 
for WR stars. Agreement is very good, except for the brief, low mass- 
loss B supergiant phase and highly uncertain high mass-loss RSG phase. 
Continuing studies into the latter are highly desirable. Typical outflow 
velocities for each evolutionary phase are also indicated in Table 4. 

Table 4 Summary of wind properties during various phases during the evolution of 
a rotating ~40M Q star (VS n i t =300 kms -1 ) at solar metallicity (following Meynet & 
Maeder 2000). Two sets of mass- loss rates are given - those adopted by Meynet & 
Maeder, M MM , and those taken from recent empirical calibrations, M omp . 



Age 


Mass 


logL 


Sp Type 


M MM 


M emp 


Voo 


Myr 


M Q 


L Q 




l(T 6 M e yr - 1 


km s" 1 


0.0 


40 


5.4 


04 V 


1.2° 


0.4 e 


3000 h 


3.7 


35 


5.5 


071 


2.l a 


2.1 e 


2100^ 


5.1 


31 


5.7 


B3I 


7.6 a 


0.7 e 


490* 


5.1 


31 


5.8 


AO I 


22. 6 


14. e 


170* 


5.1 


29 


5.8 


M0 1 


400. 6 


86/ 


3(P 


5.3 


18 


5.7 


WNL 


30. c 


26.3 


750 h 


5.3 


16 


5.6 


WNE 


25. d 


31.3 


2000 ft 


5.4 


13 


5.5 


WCE 


15. d 


28. 9 


2500^ 



(a) Lamers & Cassinelli (1996); (b) de Jager ct al. (1988); (c) Nugis et al. (1988); (d) Langcr 
(1989) corrected for clumping (Schmutz 1997); (e) Kudritzki & Puis (2000); (/) Vanbeveren 
ct al. (1998); (g) Nugis & Lamers (2000); (h) Prinja et al. (1990); (i) Lamers ct al. (1995); 
(jr) Jura & Kleinmann (1990) 
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